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Can anything from Noether's theorem 
be salvaged for discrete dynamical systems? 

Silvio Capobianco * and Tommaso Toffoli 



The dynamics of a physical system is linked to its phase-space geometry by Noether's theorem, which holds 
under standard hypotheses including continuity. Does an analogous theorem hold for discrete systems? As 
a testbed, we take the Ising spin model with both ferromagnetic and antiferromagnetic bonds. We show 
that — and why — energy not only acts as a generator of the dynamics for this family of systems, but is also 
conserved when the dynamics is time-invariant. 
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1 Introduction 

In the last three hundred years, analytical mechanics has 
turned up a wealth of profound concepts and powerful 
methods for dealing with the dynamical systems of physics. 
Can these results — of which Noether's Theorem is one of 
the best known — be applied to the dynamics of informa- 
tion systems, which today make up such a large part of our 
culture, and to what extent? 

Here one may anticipate two difficulties. On one hand, 
the above results take advantage of specific aspects of 
physics and may not be valid outside that territory. On the 
other, the methods of analytical dynamics are intimately 
tied to the calculus of the continuum, and may fail when 
applied to discrete systems such as automata and networks. 

In sum, is analytical dynamics just a "bag of tricks" for 
doing continuum physics, or does it provide templates for 
tools useful in a more general context? Does the insight 
that it offers transcend the distinction between physical and 
man-made systems (cf. [12] ) and that between discrete and 
continuous ones? 

The prototypical dynamical system that we all know and 
love, namely, physics, displays a striking property. Its ex- 
plicit evolution law has the form of a lookup table that to 
each possible state of the system — say, a 2n-tuple of real 
variables (or a vector) consisting of n positions and n mo- 
menta (where n is the number of degrees of freedom) — 
associates another 27T,-tuple — its "next state." For physical 
systems, it turns out that this table can be "compressed" 
into one that to each state associates a single real vari- 
able — the energy of that state. The first table explicitly 
expresses the dynamics as a vectorial function, while the 
second implicitly expresses the dynamics through a scalar 
function. This "compression" is nondestructive, that is, 
the full table can be reconstructed from the compressed ta- 
ble. Both the 2n-to-l encoding and the converse l-to-2ri 
encoding are achieved by simple, standardized algorithms 
(in analogy with "zipping" and "unzipping" a file). For 
this reason, the energy function, that is, the above energy 
lookup table, is said to be the generator of the dynam- 
ics, in the sense that plugging this function into a generic 
kinematic scheme (such as Hamilton's canonical equations) 
allows one to explicitly reconstruct a specific dynamics. 
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The usual way to derive the above property of energy 
relies on the fact that a physical system (both in classical 
and quantum mechanics) is imagined to evolve continuously 
in time. However, in order to have this property, being time- 
continuous is not enough — for this, the dynamics must have 
much additional structure, namely, symplectic, in classical 
mechanics, and unitary, in quantum. 

As we generalize the concept of dynamical system, in par- 
ticular to systems that, like automata of all kinds, have dis- 
crete state variables and evolve discretely in time, the ques- 
tion arises whether an energy-like quantity can be retained 
at all, and, if so, under what conditions. The first question 
we'll ask here is whether a cellular automaton can have 
something that can rightfully be called "energy," what con- 
ditions it would have to satisfy for that, and what it would 
look like. We'll then ask whether and to what extent the 
method proposed by Noether for deriving conserved quan- 
tities from the symmetries of the dynamics of a continuous 
systems may be adapted to discrete systems. 



2 Energy conservation 

We've all heard of the principle of conservation of energy: 
In an isolated physical system, no matter what transforma- 
tions may take place within it, there is a quantity, called 
energy, that remains conserved. Our belief in this principle 
is so strong that, whenever a violation of energy conser- 
vation is observed, our first reflex is not to question the 
principle, but to look more carefully into the experimental 
evidence and ask whether something may be escaping us. 
To defend the reputation of energy as a bona fide conserved 
quantity against new incriminating evidence, now and then 
physicists are even willing to stretch energy's very definition 
(see Feynman's peerless parable [S]). 

How far shall we be able to retain a viable concept of en- 
ergy as our conception of "dynamical system" keeps widen- 
ing? Can energy be conserved in a cellular automaton? 
(cf. [4 ) Conserved or not, should an energy be expected 
to be found at all in the latter? Is the idea of "energy" 
meaningful for a Turing machine or a finite automaton! In 
sum, what features does a dynamical system need to exhibit 
in order to possess a recognizable generalization of energy? 
Why should energy be conserved? 
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3 Noether's theorem 

Nearly a century ago, Emmy Noether formulated a prin- 
ciple of great generality and beauty, commonly known as 
Noether's Theorem^ [T3], which we shall state for a mo- 
ment just in headline form, ignoring the fine print: "To 
every symmetry of a dynamical system there corresponds a 
conserved quantity." According to this principle, the con- 
servation of energy (a state variable) would be the conse- 
quence of the constancy of a system's dynamical laws, that 
is, of the invariance in time of the system's structural pa- 
rameters. 

Noether's theorem plays a unifying role in our under- 
standing of mechanics; according to [13], it "has become a 
fundamental tool of modern theoretical physics." Since it 
can take advantage of many kinds of symmetries, it is very 
productive: for instance, it predicts that, if a system's laws 
are rotationally invariant, then, just because of that struc- 
tural invariance, there will be a specific state variable that 
is conserved, namely, angular momentum. 

However, when we get down to reading Noether's theo- 
rem's fine print, we discover a number of restrictive provi- 
sos. A fuller version of the contract would spell "To ev- 
ery one-parameter, continuous group of symmetries of a 
Lagrangian dynamical system there corresponds a scalar, 
real-valued conserved quantity." Thus, a discrete group of 
symmetries, such as the invariance of Conway's "game of 
life" under quarter-turn rotations of the orthogonal grid, 
won't do (bye-bye angular momentum?). Similarly, even 
though the dynamics of life is the same at every step, and 
thus time-invariant, it is not invertihle (and thus only forms 
a semigroup, not a group), and in any event the game ad- 
vances in discrete steps, instead of along a one-dimensional 
continuum (bye-bye energy?). 

4 The Ising playground 

The best way to identify the issues at play is to examine a 
specific class of models, and then map the conceptual ques- 
tion of our investigation to concrete features of a model. For 
our purposes, we shall choose a standard cellular-automaton 
realization of the Ising spin model (cf. [21 13119]). This con- 
sists of a regular array of elementary magnets, or spins, that 
can assume one of two orientations — 'up' and 'down'. 

Spins interact with one another. In the Ising model styl- 
ization of physics, long-range forces are ignored and short- 
range ones represented by the interaction of a spin with 
only its immediate neighbors: its four first neighbors, in a 
2D orthogonal array. In the simplest case, these four bonds 
(or "couplings" ) are all of a ferromagnetic nature, that is, 
they behave like rubber bands that are relaxed when the 
spins at the two ends of a band point in the same direction 
("parallel"), and stretched, when in the opposite direction 
("anti-parallel"). 

Different kinds of dynamics may be used for the Ising 
model; in order to have a unifying criterion for prescribing 
a dynamics, one starts with assigning one unit of a notional 
"potential energy" to stressed bonds, and zero to unstressed 
ones. The class of dynamics we choose here, the so-called 



microcanonical Ising model, are strictly deterministic and 
invertible: on a given step, a spin will flip (that is, reverse its 
orientation) if and only if doing so will leave the sum of the 
potential energies of the four surrounding bonds unchanged. 
Note that, in the ferromagnetic case, this happens when 
two of the neighbors are parallel to the given spin and the 
other two anti-parallel. In this case, in fact, we start with 
two stressed bonds. After the flipping of a given spin, its 
parallel neighbors will have turned anti-parallel and vice 
versa, so that the two stretched bonds relax at the same 
time as the two relaxed ones stretch, and consequently the 
overall potential energy stored in these bonds is two units 
both before and after the step. In all other cases, since no 
flip occurs, potential energy is conserved trivially. 

A final stipulation. The above rule would become incon- 
sistent if one attempted to update two adjacent spins at 
the same time. In fact, if both spins are instructed to flip, 
the assumption on the part of either spin that the shared 
bond's stress status will thereupon be complemented fails 
(as its stress status will remain unchanged), and energy is 
no longer conserved. A standard solution is to treat the ar- 
ray of spins as a checkerboard consisting of two intermeshed 
subarrays called even and odd (the red and black squares 
of the board) according to the parity of x -\- y, where x 
and y are the spatial coordinate of a site. A system state 
s consists of an ordered pair {qA,qB) of configurations of 
different parity (i.e., one even and one odd, in either or- 
der). Together, an even and an odd configuration fill up 
the whole array; however, their order in the pair is relevant, 
that is, state {qA,<lB) is distinct from {qB,qA)- Depending 
on the context, it may be convenient to call the two con- 
figurations of a pair either the "past" and the "present" 
or, alternatively, the "present" and the "future" (the latter 
abbreviated, when subscript, as "pres" and "futr"). 

The Ising dynamics is a second-order recurrence relation 
of the form 

1x,y ® fiQx-l,yjQx+l,yTlx,y-lTlx,y+l) — Qxly ■ (1) 

This is a conditional permutation: a spin q^.y in the 
past (time t — 1) is mapped into a spin q^^y in the fu- 
ture (time i + 1) by xORing it with an enable bit — 
the output of the binary function / when given as argu- 
ments the present values (time t) of the first-neighbor spins 
qx-i,y,qx+i,y,qx,y-i,qx,y+i, which bclong to the comple- 
mentary subarray and thus affect the updating of qx,y but 
are not affected it. Depending on the value of the en- 
able bit, this updating either leaves a spin unchanged or 
complements it — both invertible operations. The above 
"shifting-frame" mechanism — of the form ((/past, 9pros) '-^ 
(gpios, 9futr), where the right item qp^os of the old state shifts 
to the left unchanged, pushing out the left item gpast, and 
a new item (/futr, computed according to (|T]), shifts in from 
the right — is called in numerical analysis leapfrog updating. 

The macroscopic behavior of this kind of model is trivial 
in one dimension; conceptually productive in two dimen- 
sions, as it provides by very simple means an insight into 
the physics of phase transitions; and increasingly rich and 
challenging in more dimensions. 
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5 Why call it energy? 

Setting physical motivation aside, we have in the Ising spin 
model an abstract symbolic dynamics consisting of a two- 
dimensional array over the binary alphabet {0, 1} (for what 
was 'down' and 'up'; it will still be convenient to call "spin" 
the state variable at each site) . We distinguish between even 
and odd sites according to the parity of the sum x+y of their 
coordinates. We update even and odd sites on alternating 
steps according to the rule (stated the rule in a way that 
obviously generalizes to any number of dimensions) that 
"A spin u at (x, y) flips (i.e., complements) if and only 
if its first neighbors (those at {xztl,y) and (x, y ± 1)) 
are evenly divided between state u and its complement 
u." 

Note that we had no use for an "energy." However, if as an 
afterthought we choose to call energy a distributed quan- 
tity whose value, for each pair of adjacent first-neighbor 
spins, is if these two spins are equal and 1 if different, 
then it follows that energy is conserved by the above dy- 
namics. Also note that this energy is but the length of 
the boundary between 'up' and 'down' domains. While the 
magnetization — number of spins up minus number of spins 
down — may change with time, that length, and thus the 
energy, remains constant. Is this use of the term "energy" 
just physics nostalgia — like calling a British colony "New 
Hampshire"? Or is there something more to it? 

For most people, energy is a fungible resource of the form 
E = T + U , that is, stored either in the inertia of a moving 
object ("kinetic energy" T) or in the state of compression 
of an elastic medium such as a spring ("potential energy" 
U) . Such a concrete presentation of energy will do for many 
practical applications. 

But if you look at what energy is really supposed to be in 
physics, you'll discover an abstract quantity that is totally 
unconcerned with the nature of the objects and materials 
that make up a dynamical system, and instead totally ab- 
sorbed with the topological structure of the system's trajec- 
tories. In fact, the total energy of a system may be defined 
as 

1. A real-valued function of the system's state, 

2. that is additive, (2) 

3. and is a generator of the dynamics. 

Incidentally, there is no mention of "energy conservation" 
in all this. 

Note. We wish we could give a precise reference for the above 
definition. In fact, here we rely on disparate loose hints culled 
and integrated with some effort from the literature. Physicists 
are supposed to know what energy is, and they don't need — 
perhaps have good reasons not to want — a final definition (of. 
the Feynman parable mentioned in So, not only college text- 
books, but even such beloved conceptually-minded references as 
Arnold[T], Goldstein[6], and Lanczos[7] proceed by giving exam- 
ples of energy in increasingly elaborate contexts, and then bring- 
ing one's attention to increasingly more abstract properties of 
energy. Definitions — hardly ever. 

Let's verify that the "energy" of the Ising model agrees 
with our abstract definition. In order to do so, we'll be 



forced to re-examine the definition and pore over the fine 
print that by tacit agreement stands behind such terms as 
"state," "additive," and "generator." 

5.1 Additivity 

As we've seen, what we proposed to call Ising "energy" of 
a state is simply the number of streteched bonds in that 
state. Incidentally, since the Ising model consist of an in- 
finite array of sites, the state of the system consists of an 
infinite array of spins, and virtually any state will have in- 
finite energy. For a function to be additive, here the physi- 
cist means that (a) it is meaningful to partition the system 
into separate subsystems such that each of these subsytems 
has its own well-defined state, and the global state is the 
composition, or Cartesian product, of the states of these 
subsystems; (b) the function is defined on the state of each 
of these subsystems, and (c) its value for the whole system 
is the sum of its values for the individual subsystems. 

If a system could be partitioned into subsystems or 
"blocks" that are totally independent of one another, each 
one having a well-defined energy independent of the state 
of the other blocks, additivity would be trivial to achieve — 
just define the energy of the whole system as the sum of 
the energy of its subsystems — but by the same token would 
be a vacuous concept. In general, however, there will be 
some residual coupling between blocks, and even if one knew 
the precise state of a block one would be able to assign an 
energy to it only to within an upper bound and a lower 
bound, which refiect the lack of knowledge of the block's 
environment's state. Fortunately, in many cases, when in- 
teractions are local and the blocks are appropriately chosen, 
this uncertainty is roughly proportional to the length of a 
block's boundary — where it interfaces with its environment, 
and thus the relative uncertainty is bounded above by the 
boundary /volume ratio for a block, decreases with block 
size, and vanishes in the limit of arbitrarily large blocks. 
This is what "additivity" means to the physicist, and it is 
consistent with the approach taken by topological dynam- 
ics. 

In the Ising model, sites interact with one another 
through an unbroken mesh of couplings; there are no strictly 
independent subsystems. However, the stipulation made at 
the end of [JH whereby the odd mesh is kept constant when 
the even mesh is updated (and vice versa) — and thus ef- 
fectively acts as a parameter rather than a state- variable — 
enforces a discipline by which energy variations correspond- 
ing to local variations of a configuration can be determined 
exactly. We shall see this at the end of i )5.2l and use it to 
full advantage in ^Zl 

5.2 Generator of the dynamics 

A generator of the dynamics is a function of the system's 
state through the knowledge of which one can reconstruct 
the system's dynamics in an explicit, vectorial form (cf. 
^ — of course only up to an isomorphism. It is under- 
stood that this reconstruction should be achievable by an 
algorithm given once-and-for-all, independent of the specific 
system (otherwise, one could "cheat" and hide information 
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about the system in the algorithm itself — big deal!). Sup- 
pose the function is a "Hamiltonian" that for each state 
{<1,P) = {qi,---,qn,Pi,---,Pn) yields its energy H{q,p). 
Our goal is to determine where any state will go. 

For a continuous dynamics as we imagine in the present 
example one would want the instantaneous "direction and 
rate" at which the state progresses, that is, the 2n deriva- 
tives {q = — III other words, we want to arrive at 
vectorial mapping, of the form 



{q,p) ^ {q,p), 



(3) 



even though mere evaluation of the function at {q,p) will 
perform the mapping {q,p) i— >■ H{q,p) and thus only give 
us a a scalar. However, if we performed repeated samplings 
of H in the vicinity of {q,p) (after all, we have the entire 
function H and we can apply it as many time as we want) 
so as to get a "sense of direction," or, equivalently, in the 
continuum case, if we could get the derivatives of H at that 
point, we would get more information associated with that 
point than its mere energy H{q,p). In Hamiltonian me- 
chanics, mapping ([3]) is provided by the standard canonical 
equations 

. dH . dH 



q = 



dp 



P = 



dq 



(4) 



Note that, since q and p are actually abbreviations for n- 
tuples, here we have 2n equations yielding a vector consist- 
ing of n pairs {qi,Pi), (i = l,...,n). 

It will presently become obvious that in the Ising model 
the energy function does generate the dynamics — and by 
means of a very simple recipe. In a discrete Hamiltonian 
dynamics, a state is no longer a "position/momentum" pair 
{q,p) as in the continuous case, but an ordered pair of con- 
figurations {qo,qi), as explained at the end of The rules 
of the game are that one may propose any such state, as 
many times as one wislies, and each time get as an answer 
the energy i?((?past, ^pros) for that state. Isn't this number, 
in the range {0, 1,2,..., oo} (and almost always oo) much 
too little to determine the infinite qfutr configuration? The 
answer lies in the principle of virtual displacements. That 
is, given (/past and gprcs, one may propose any conceivable 
candidate for qfatr, ask for the energy of the resulting new 
state ((7pics, 9tuti)7 compare it with the old energy, and on 
the basis of this comparison decide whether to accept or 
reject the candidate. Just as in the "hot, hot, cold, cold" 
game, in the game of "battleships" — or, for that matter, in 
biological evolution — one is never explicitly shown the tar- 
get. Instead, one must describe an object, and will only be 
told ("yes," "no") whether that is target, or possibly given 
an intermediate scalar value ("hot", "warm", etc.) that 
indicates how far one's proposal is from the target. In prin- 
ciple, one can arrive at the solution by going through all 
possible candidates. The point of the game is to speed up 
the search by dynamically arranging the proposals so that 
the feedback will give one a "sense of direction." 

In the Ising case, we want to find a qtutr such that the en- 
ergy difference AE = iJ((7prcs, gfutr) - -ff (^past, gprcs) is zero 
not only globally but also locally, that is, at any scale down 
to the radius of the first-neighbor neighborhood (the range 



of interspin "interactions" ) . Forget about getting finite en- 
ergies by taking finite configurations by blocking (i i|5.1|) : just 
take the 32 infinite configurations of the form 





w 





(5) 



with c denoting a "center cell" in the past and n,s,w,e 
denoting its four neighbors in the present. This infinite 
state will have a finite energy ranging from through 4. 
Let's make a table of energy for the 32 possible collective 
values of those five sites, 



c nswe 



0000 
0001 
0010 
0011 
0100 
0101 
0110 
0111 
1000 
1001 
1010 
1011 
1100 
1101 
1110 
1111 



c nswe 





1 

1 

2 1 

1 

2 1 

2 1 

3 

1 

2 1 

2 1 

3 

2 1 

3 

3 

4 



1 0000 
1 0001 
1 0010 
1 0011 
1 0100 
1 0101 
1 0110 
1 0111 
1 1000 
1 1001 
1 1010 
1 1011 
1 1100 
1 1101 
1 1110 

1 1111 



Hf 



4 
3 
3 

2 1 

3 
2 1 

2 1 

1 

3 

2 1 
2 1 

1 

2 1 
1 
1 




(6) 



where H is the "Hamiltonian energy" of the Ising model for 
the entire global state ([5]), / the "conditioning" function 
that appears in ([1]), and c', as we shall see in a moment, the 
effective future value of the center site. As we recall from 
iJH the Ising dynamics is "flip if and only if that leaves the 
energy of the neighborhood unchanged." Therefore, for any 
state of the form ([5|) , we may propose as candidates for the 
new state either (a) the same state, i.e., with c' = c, or (b) 
the one where c has flipped, i.e., where c' — c. By using 
table ([6]) two times, i.e., for the two c' candidates c and 
c, we see that the energy is (obviously) unchanged in the 
first case, and may have changed in the second case. If the 
energy is the same in both cases, then the actual value for 
c' (that is, the one prescribed by the Ising dynamics) will 
be the default value c; otherwise, it will be c. This is the 
value given in column c' of (jG]). 

Finally, we must show that this new value c', given the 

n 

old value of the entire neighborhood "'ge, is independent of 
all the others sites of the global state — those assigned in 
global state©. In fact, by the definition, in the Ising model, 
of energy as a sum over bonds (or adjacent-spin pairs), the 



energy spanning, say, the "extended patch" 



is the sum 



of the energies spanning the two elementary patches wee 
and nar-^ in other words, energy is strictly additive, and the 

local dynamics table (|n|) can be applied separately to every 
elementary patch. 

Thus, the "energy" of the Ising model satisfies the ab- 
stract, "model- independent" definition of energy given in 
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(g]). ... Or does it? 

6 Raising the kinematics bar 

Indeed the Ising energy (a) is a real-valued function of the 
system's state; (b) it is additive, since the alternating up- 
dating scheme guarantees that the accounting of energy can 
be done independently on every elementary patch; and (c) 
a simple rule allows us to reconstruct the dynamics by sam- 
pling the energy function in the (topological) neighborhood 
of a state. 

However, point (c) is in this case vacuous, since here we 
are dealing with a single system — the Ising model — and 
there is no indication of how one would go about gener- 
ating the dynamics of other systems using other energy 
functions. Consider a telegraph channel with a single in- 
put key, pressing which will invariably transmit the text of 
the Kamasutra. In spite of the length and the intrinsic in- 
terest of this text, the capacity of the channel is obviously 
log 1=0 bits. In a similar way, the information contents of 
-falsing as a function, for the sake of specifying a dynamics, is 
zero, because it is a single "key" : the Ising dynamics is the 
only dynamics we have in mind, it is known, and we might 
as well have stored a copy of it at the channel's receiving 
end. What one needs is a generic function scheme — a look- 
up table "blank form," as it were — which, when filled with 
specific contents and processed by an algorithm given once 
and for all, will allow one to describe any one of an agreed 
upon class of dynamical systems. 

This issue is addressed by physicists by distinguishing be- 
tween "dynamics" and "kinematics," where the latter term 
is used in a somewhat specialized sense. That is, a kine- 
matics is a framework — a set of rules — that characterizes a 
whole class of dynamical systems of a similar type or "tex- 
ture," while a dynamics is any one of the systems sharing 
that kinematics. Until we place the Ising model within a 
larger class of models, and show the existence of an energy- 
function- valued variable — whose values will be specific en- 
ergy functions each of which will in turn allow a fixed algo- 
rithm to generate the different dynamics of that class — until 
then we will not have a bona fide generator of the dynamics. 

There is a natural and a very simple generalization of 
the Ising model (used, specifically, in the modeling of spin 
glasses) that yields a class of sixteen dynamics. We'll show 
that for that class the energy function plays a nonvacuous 
role as a generator of the dynamics. 

In fJH we called ferromagnetic a bond between adjacent 
spins that is unstressed when the two spins are parallel and 
stressed when antiparallel ("type 0" bond). We now in- 
troduce a bond of "type 1," or anti ferromagnetic, which is 
stressed when the spins are parallel, and unstressed oth- 
erwise. With bonds of two kinds thus available, and four 
bonds surrounding a spin, there are 2^* = 16 possible struc- 
tural environments by which a spin can find itself sur- 
rounded. We shall retain our definition of Ising energy as 
the sum of the energy of the bonds, recalling only that the 
stress u -f-^ ti of a bond of type b between a pair of parallel 
spins u and v will u © 6 © w, rather than just it © u as in the 



original Ising model. As part of our kinematics — not as a 
peculiarity of an individual dynamics — we shall retain the 
Ising prescription that a spin will flip if and only if such a 
move leaves the energy unchanged for a given assignment 
of the four surrounding bond types. However, within this 
kinematics, we shall for the moment restrict our attention 
to systems whose structural parameters (those determining 
a dynamics, and in our case consisting of the assignment of 
bond types) be time- and space-invariant, that is, for each 
specific dynamics of our kinematics, the behavior of a spin 
shall not depend on where the spin is located and on what 
time it is, but only on the values of its neighbors. 

With the latter restriction, we are left with as many 
global assignments of structural parameters as there are 
choices of local assignments — that is, the 16 distinct as- 
signments of four bond types around a spin. For each of 
these structures we must tabulate, by running over all pos- 
sible states, the corresponding dynamics (a function from 
old global state to new global state) and the corresponding 
energy function (a function from global state to an inte- 
ger). Observe that two distinct structures may well end 
up yielding the same dynamics or the same energy func- 
tion. Our goal will be achieved if (a) no two dynamics are 
specified by the same energy function (in other words, the 
energy function distinguishes between dynamics), and (b) 
every dynamics has at least one energy function that spec- 
ifies it. Since, by construction, energy is additive down to 
the scale of individual four-bond "patches," the local "vari- 
ational principle" employed in H5.2I will then allow us not 
only to indirectly denote, but to explicitly generate the cor- 
responding vectorial-form (that is, state-to-state) dynam- 
ics. 

This is indeed the case. Since we are presenting here an 
original approach, it is important not to deceive ourselves on 
delicate aspects of an argument. For this reason, as a "san- 
ity check" we have explicitly tabulated dynamics and energy 
function for all 16 elements of our generalized Ising kinemat- 
ics on a finite, toroidally wrapped-around model. (We used 
just a 4x4 torus; though small, this size is large enough to 
rule out spurious degeneracies. Full scripts, output data, 
and statistics are available at www. ioc . ee/~silvio/nrg/',) 
A formal proof is subsumed within the properties of the 
numerical integration scheme of Fig.[Tl which constitute 
a proof in the more general context of space- and time- 
dependent dynamics. In fact, in the next section we shall 
remove the constraint that the structure — and the atten- 
dant dynamics — be spacetimc uniform. 

7 Space- and time-dependent dynamics 

Energy in these generalized Ising models is a function of the 
global state — a function ultimately specified by the struc- 
ture parameters of the specific model. If the parameters 
are allowed to change from one step to the next, then 
this function is in turn a function of time. In fact, the 
type of the bond joining two given spins will depend on 
"what time it is" at the moment of performing a transi- 
tion {qA,qB) {lB,<lc)- The "new state" ((7b,9c) at the 
end of this transition reappears as the "old state" at the 
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beginning of the following transition, ((7b,<Zc) ^ (qCiQu), 
but now it is a difFerent time, and the energy of the same 
state (of spins) has to be re-evaluated accordingly to the 
new assignment of structural parameters. 

This arrangement of things, which is essentially a numer- 
ical integration scheme, is graphically illustrated by Fig.[T] 
There, the integers along the time axis represent discrete 
time instants, those along x or y, discrete spatial positions. 
A spin is represented by a thick vertical line, interrupted 
by a conditional-permutation gate (©) every time that it is 
the spin's turn to be updated. A system state is an ordered 
pair of configurations associated with two consecutive in- 
teger times. This state may conveniently be labeled with 
the "half-time" in between; that is, one may denote by ^ 
the state consisting of the configurations at and 1, and so 
forth. Steps of the dynamics are nominally associated with 
integer values of time. So, for instance, state l| is the "new 
state" for step 1 but the "old state" for step 2. Horizontal 
arcs of unit length (thin lines) represent interspin bonds as 
at the time of a step; letters may be used to speciiy the na- 
ture of a bond; so r denotes the nature of the bond between 
spins 1 and 2 at step 1, while s denotes the nature of the 
same bond at step 2. 




Figure 1: (a) Spacetime spin grid {y not shown), (b) Detail of 
full grid, with both x and y. 

The value of a bond remains steady during spin updating 
(at integer times) but may change between updatings (at 
half- integer times). On this occasion, the energy of a state 
is accordingly re-evaluated. Thus, even if energy is con- 
served during spin updating, in a time-dependent system 
the energy may still change when the bonds are altered. 

Thus, it will generally happen that energy changes with 
time. With the present full class of dynamics (i.e., one 
which allows arbitrary spacetime variation of the structural 
parameters) of the given generalized Ising kinematics, we 
arrive at the conclusion that energy is necessarily conserved 
if the spacetime structure is time-invariant (because the 
transitions themselves are by construction energy conserv- 
ing, and with time-invariant bond types no re-evaluation 
of the energy need occur). On the other hand, when the 
spacetime structure is not time invariant, at least some of 
the dynamics are not energy-conserving. We can thus right- 
fully say that, in the present context, "Ising energy" as a 
whole (not the energy of a specific Ising dynamics) is the 
quantity that is conserved because of the time-invariance 
of the dynamics. 

For this class of dynamics, then, in spite of the dis- 
creteness of the dynamics, a non-frivolous transliteration of 
Noether's theorem applies with full force, yielding energy 
conservation from time invariance. 



8 Conclusions 

We've shown that certain aspects of Noether's theorem ap- 
ply to dynamical systems beyond ordinary analytical me- 
chanics, and why they do. We were able to shed the require- 
ment that the dynamics be continuous; still, some form of 
second-order discipline seems to remain essential for a sym- 
metry to give rise to a conservation law. This is just the 
beginning of what promises to be a productive line of re- 
search. 
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